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Abstract. 

It is shown that Lame's equation -j^X + K^cn^(z, — can be reduced to 

hyper-geometric equation. The characteristic exponents of this equation are expressed 
in terms of elementary functions of the parameter /t. An analytical condition for 
parametric amplification is obtained. 
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1. Introduction 

Lame's equation 

—X+{Ri + R2sn\z,K))X = Q, 1 

where sn{z, K) is the standard Jacobian elliptic sine function of modulus K ||, and 
i?2 are numerical parameters, is common in many branches of mathematical physics. 

Attention to this equation has been called recently by the fact that it plays an 
important role in certain problems of particle physics and cosmology of the Very Early 
Universe. In particular, it was shown [1] that in the theory of two fields (f) and x with 
potential V{(f), x) = ^ + ^-^^ the equation of motion of the field x iii the Heisenberg 
representation can be reduced to Lame's equation under certain assumptions. More 
concretely, for the Minkowski space-time in the linear approximation for the field x, it 
is found that after expansion of this field over the eigenf unctions of the Laplace operator, 
the equation of motion for an eigenfunction corresponding to a particular wavenumber 
has the form 

where cn{z^ ■^) is the Jacobian elliptic cosine function of modulus K = the 
parameter k is proportional to the wavenumber, and the "time" z is proportional to 
the ordinary Minkowski time coordinate. In fact, the same equation is valid in the 
expanding Universe in a certain regime provided that the eigenfunctions are rescaled in 
a proper way, and the coordinate z plays the role of so-called conformal time. 

Along the real axis, cn{z, is a periodic function with period T = 4K(-^) = 

^ ^'^^ ^ 7.416, where K(i^') is the complete elliptic integral of modulus K of first type 
and r(x) is the gamma function. It follows from general theorem that the equation (2) 
must contain solutions in the form Xi^2 = e'^^''^^P{z), where P{z) is a periodic function 
of the period T, and the coefficients /ii_2 are called characteristic exponents. If one of 
these coefficients is real and positive, the corresponding solution describes an exponential 
growth of the amplitude of the eigenfunction X. In modern theories of matter creation 
in the Universe [2-4] this growth is interpreted as a production of "particles" of the field 
X- The rate of production of the "particles" is determined by the values of characteristic 
exponents, and therefore the calculation of these exponents is very important for such 
theories. Usually the calculation of characteristic exponents is performed by numerical 
means, and only a few cases are known with analytical solutions. In particular, the 
characteristic exponents have been calculated analytically in the paper [1] for the case 
g_ _ n(n+i) integer). 

In this note we would like to point out that the special case of the equation (2) 
with A; = 0: 

i^X + «W(..-L)X = 0, 3 

I The elliptic sine function sn{u,K) and the elliptic cosine function cn(u,K) are defined with help of 
the elliptic integral u = , , Then, sn(u,K) = sin((/)) and cn{u,K) = cos(0). 
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(k = ^) can be reduced to the hyper-geometric equation. This allows for exact 
calculation of the characteristic exponents for that important case ||, representing them 
in a remarkably simple form (see eqs. 21, 25 below). 



2. Reduction of the equation to the hyper-geometric equation and the 
characteristic exponents 

Let us consider the equation of form (3), and make the following change of the 
independent variable: 

y = cn^(^, _). 4 
Then, we use the well known relation 

^cn{z, K) = ^(1 - cn^{z, K)){K'^ + K^cn^{z, K)), 5 



where i^' = y (1 — K'^). Taking into account that in our case we obtain 

where the + sign (— sign) should be taken if y increases (decreases) with z. Also, we 
have 

Substituting the last operator into equation (3), we see that the factor y^/^ cancels and 
that this equation is transformed to the standard hyper-geometric form: 

yil - y)£^X + (7 - (a + + l)y)-^X - a(5X = 0, 8 



where 7 = |, and a,l3 = |(1 ± y (1 + 8^^)). The general solution to eq. (8) can be 
written in terms of two elementary solutions valid in the vicinity of the singular point 
2/ = 0: 

X = CiCpi + C202, 9 

where 

is the Gaussian hyper-geometric function, and 

h = y'-^'F{a + 1 - 7, /3 + 1 - 7, 2 - 7; ?/) • 11 

Alternatively, the solution to eq (8) can be written in terms of two elementary solutions 
valid in the vicinity of the singular point y = 1: 

X = C303 + C404, 12 
§ Which corresponds to a long wave approximation from viewpoint of particle physics and cosmology. 
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where 

03 = F(a,/?,a + /3 + l-7;l-y), 13 

and 

04 = (1 - y)^-"-^i^(7 - «,7 - /3,7 + 1 - a - /3; 1 - y). 14 
The pairs of solutions 01, 02 and 03, 04 are connected by the well known relations (e. 
g- [5]) 

01 = A03 + S04, 02 = C03 + D04, 15 

r(7)r(7-5) ^ _ r(7)r(^ - 7) 



r(7-a)r(7-/3)' r(a)r(/3) 

r(2 - 7)r(7 - 5) r(2 - 7)r((^ - 7) 

r(i-a)r(i-/3)' r(a + i-7)r(/? + i-7)' 

where r(a:) is the Gamma function, and S — a + p. 

Let us point out that the transformation (4) is singular at the points where the 
function y is equal to zero or unity — ^ when y ^ 0,1). A simple analysis shows 
that the coefficient C2 must change its sign at the points Zi defined by y{zi) = 0, and the 
coefficient C4 must change its sign at the points Zj satisfying y{zj) = 1. To show that, 
let us consider the behavior of the function X near the points Zj. Near these points, we 
can approximately write 

^ ~ ci + C2y'/^ 18 
and taking into account equation (6), we have 

— p=i -, 19 

dz y/2 

provided ^ < 0, and 

dz V2' 

provided ^ > 0. The function X and its derivative with respect to "time" z must be 
continuous functions of z. Therefore, the coefficient C2 must change its sign at the points 
Zi. Similar arguments can be used to prove that the coefficient C4 must change its sign 
at the points Zj where y{zj) = 1. 

Thus z changed over half a period of cn{z, ■^), a new decomposition of the solution 
of eq. (8) should be made: 

X = Ci0i + C202, 21 

where in general the coefficients £1^2 do not coincide with the coefficients Ci^2- The rules 
for changing of these coefficients follow directly from the arguments mentioned above if 
one uses eqs (9-15) taking into account the explicit form of the connection coefficients 
(16,17). It is straightforward to obtain the relation: 



? = t)c^, 



22 
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where the components of the matrix t* have the following explicit form: 

tl = tl = V2cos7r{a- P), 23 

' r(i - a)r(i - /3)r(i + « - 7)r(i + /? - 7) ' 
' r(a)r(/?)r(7-a)r(7-/?)' 

The eigenvalues of the matrix are 



Ai,2 = "\/2cos7r(a — /5) ± Jcos2TT{a — (3), 26 



and a - /3 - 
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are ^ 2 

— In pi_2 = — ln(-\/2 cos7r(a — /?) ± -J cos 27r(a — 27 



Obviously, the multiplicators pi,2 are equal to 25 and the characteristic exponents 

/^l,2 ^— rx,^ y 

In actual applications it is very important to know under which conditions a 
particular solution to eq. (3) experiences parametric amplification, that is under 
which conditions its amplitude increases when the "time" z changes over the period of 
cn{z, -^). To characterize the parametric amplification, we introduce the real quantity 

Ii{k,) = Max{Re{ni^2)), 28 

Obviously, > is the condition for parametric amplification. It is easy to see that 
it is satisfied when Ai,2 are real, that is when 

n--<a-d<n + -, 29 
4 4 

where n is an integer > 1. In terms of the parameter k, these inequalities can be 
rewritten as 

^n{2n - 1) < K< ^n{2n+ 1). 30 

The expression for the quantity ft follows from eq. (27) and the definitions of the 
parameters a and f3 and the expression for the period T = 4K(-^) = 



2^ , r ATI + 8^2 / TTv/TTS 
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m = ^ In m cos + ^cos c^'^ n 31 

The quantity attains its maximal value at 



where its value is 



'2n2--, 32 



/imax = In (1 + v^) ^ 0.2377. 33 



.4/ 

Note that the same value has been obtained in [1] in the asymptotic limit k — >• cxd. 
II To obtain the eigenvalues Ai,2, we use the weU known relations r(a; + 1) = xY(x) and r(l — x)V{x) 

TT 

sin {t^x) 
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3. Discussion 



We were not able to find the simple formulae derived in the standard reference books 
on the Lame's equation [6-8]. However, a very similar transformation between another 
equation of Lame's type and a hyper-geometric equation has been discussed recently by 
Clarkson and Olver [9] (see also [10]). They show that the hyper-geometric equation 

and Lame's equation 

— * + SGaWiu)"^ = 0, 35 

where W{u) is the Weierstrass elliptic function with parameters g2 — 0,93— 33/2 ^g , are 
related to each other by the transformation 



V3+i-{i-ty 



cn{u, K) = ^^—^ 36 



where K — \J\-\- Obviously, the characteristic exponents of equation (29) can 
be obtained by a method similar to that described above. In general, it would be very 
interesting to find a general solution to the following problem: under what condition can 
Lame's equation of the general form be transformed to a hyper-geometric equation? The 
solution could be applied in many problems of modern particle physics and cosmology. 

Finally, it is interesting to note that in principle the same line of argument could 
be apphed to Lame's equation of the general form. It is well known that after the 
change of variable yi = cn?{z,K) Lame's equation is reduced to a particular form of 
the Heun's equation. Then the calculation of the characteristic exponents is reduced to 
solving of the connection problem for Heun's equation between the elementary solutions 
corresponding to the singular points y\ — and yi = 1 (e.g. [11]. This problem can 
indeed be solved [11], but the connection coefficients analogous to coefficients (16,17) 
are now expressed in terms of complicated series which look rather difficult for analytic 
treatment of the general case. 



acknowledgments 

I am grateful to A. Illarionov, A. Kirillov and A. Starobinsky for remarks. This work 
has been supported in part by RFBR grant N 00-02-16135 and in part by the Danish 
Research Foundation through its establishment of the Theoretical Astrophysics Center. 



References 

1 Green P, Kofman L, Linde A and Starobinsky A 1997 Phys. Rev. D 56 6175-6192 



On Lame 's equation of a particular kind 



7 



2 Dolgov A and Kirilova D 1990 Sov. Nucl. Phys. 51 273; Traschen J and Brandcnborgcr R 1990 

Phys. Rev. D 42 2491; Shtanov Y Traschen J and Brandenberger R 1995 Phys. Rev. D 51 5438 

3 Kofman L Linde A and Starobinsky A 1994 Phys. Rev. Letters 73 3195 

4 Kofman L Linde A and Starobinsky A 1997 Phys. Rev. D 56 3258 

5 Caratheodory C 1960 Theory of Functions v 2 (New York: Chelsea Publishing Company) 

6 Whittaker E and Watson G 1996 A Course of Modern Analysis (Cambridge: Cambridge University 

Press) 

7 Batcman H and Erdelyi A 1955 Higher Transcendental Functions v 3 (McGraw-Hill Book Company) 

8 Arscott F 1964 Periodic Differential Equations (Pergamon Press) 

9 Clarkson P and Olver P 1996 J. of Differential Equations 124 225 

10 Kamke E 1975 Differentialgleichungen Losungsmethoden and Losungen v 1 (Chelsea, New York) 

11 Schafke R and Schmidt D 1980 SIAM J. Math. Anal. 11 848 



